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ABSTRACT 

High-energy photons emitted from gamma-ray bursts (GRBs) are subject to pair-production inter- 
actions with lower energy photons, leading to an effective optical depth. In this paper, we estimate 
the opacity resulting from photon fields located at various distances from long GRB sites: that of 
the binary companion to the massive stellar progenitor, that of the star- forming molecular cloud 
containing the GRB, and the total photon field of the host galaxy. The first two photon fields are 
found to be transparent for most reasonable sets of assumptions about these systems. In the case of 
galactic radiation fields, we have performed several numerical simulations to calculate the expected 
opacities for 2 different line-of-sight geometries through the host galaxy, and include a full accounting 
of the infrared radiation produced by the absorption and re-radiation of starlight by dust. The optical 
depth for GeV gamma-rays, due to direct starlight is less than unity for all host galaxies. At higher 
energies, >10 TeV, a spectral cutoff can occur due to the rapidly increasing number of mid- to far-IR 
intra-galactic photons reradiated by dust. Photons in the extragalactic background light therefore 
remain the only relevant source of photon-photon opacity for ongoing GRB observations with Fermi 
LAT, and potential future detections with ground-based gamma-ray telescopes. 
Subject headings: gamma-ray burst: general 



1. INTRODUCTION 

Recent observations by the Fermi satellite have con- 
firmed that both long and short gamma-ray bursts 
are capable of emission at energies > 10 GeV. 
The rest-frame energies inferred for gamma rays from 
long-dura tion gamma-ray bu rst (GRBs) 080916C and 
090902B (lAbdo et al.l [2009allbl) and short GRB 090510 
(|Omodei et al.l |2009[ ) are all in excess of ~60 GeV. 
These findings raise the question of whether emis- 
sion might occur at still higher energies, and these 
gamma rays might be observable by Fermi, or by 
ground-based instruments such as MAGIC, VERI- 
TAS, and HESS. Electron-positron pair-production 
from int eractions b etween gamm a rays and UV-IR 
photons (iGould fc S chreder 1 9671: iStecker et al.l 119921: 
IMadau fc Phinnevlll996t IGilmore et al.ll2009[) is an im- 
portant effect at these energy scales, modifying the spec- 
tra of distant extragalactic sources and shrouding the 
highest-energy emission. This effect could thus be a 
powerful techniqu e for understanding cosmological back- 
ground radiation ( IGilmore et aLll2010l: IAbdoll2Q10[ ) 

The photon-photon optical depth per unit path length 
for a gamma ray of energy E is 
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where n{E') is the local density of target photons 
at energy i?', and 6 is the angle of interaction. 
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Here a(E,E',d ) is the cross section for interaction 
(jMadau fc Phin ncv 1996), which peaks at twice the min- 
imum energy for pair creation 
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At energies below 1 TeV, it is the optical, UV, and 
near IR photons that are responsible for attenuation. 
Higher energy gamma rays have sufficient energy to pro- 
duce electron positron pairs in interactions with longer- 
wavelength IR photons. 

In general, attenuation at a cosmological scale is ex- 
pected to dominate the pair-production opacity, as the 
strength of the effect scales with path length through 
the intervening photon field. Analyses of the interstel- 
lar radiation fields within the Milky Way find signifi- 
cant attenuation of high energy photons only at sev- 
eral tens of TeV, with negligibl e absorption at lower 
energies (iMoskalenko et al.l 120061 ) . There is strong ev- 
idence, however, that long-duration GRBs are associ- 
ated with deaths of massive stars (jWooslev fc Blooiiil 
120061) ■ and therefore occur in star-forming regions of 
galaxies, and that the hosts of these objects often 
exhibit considerable star-formation rates ('Bloom et al.l 
[2002; Christcnscn et al. 2004; Fruchtcr ct al. 2006), and 
therefore significant number densities of UV and opti- 
cal photons. The prevailing theory for these bursts is 
that they are associated with some subset of type Ic su- 
p ernova e resulti ng from the deaths of short-lived stars 
CWoosle v 1 99 31: llwamoto et"all 119981 : iMacFadven et al.l 
I200H) . 

In this paper, we examine the optical depth for gamma- 
rays from a GRB produced by interactions with "lo- 
cal" photons within a host galaxy. Loosely speaking, 
the opacity from pair-production is proportional to the 
product of path length and density of target photons at 
the relevant energy scale. If the line-of-sight to the GRB 
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passes through a region of intense flux from nearby stars, 
the resulting optical depth could conceivably be equal to 
or greater than that from intergalactic background pho- 
tons. Attenuation effects from photon fields on several 
different length scales could be pertinent. In particular, 
we will discuss here absorption at sub-parsec scales from 
the photon field of a binary companion to the GRB pro- 
genitor, parsec-scale absorption within the star-forming 
region of the progenitor site, and finally opacity due to 
the kiloparsec-scale galactic disk of the host galaxy. In 
^ we present some simple estimates of the opacity aris- 
ing on each of these scales. In [J3l we describe a numeri- 
cal calculation performed to calculate the attenuation for 
star-forming galactic disks, followed by our conclusions. 

2. GAMMA-RAY ATTENUATION AT GALACTIC AND 
SUB-GALACTIC SCALES 

As a simple preliminary approximation, an estimate of 
r-y can be found for a given region by the product 

« nuvRa, (3) 

where R is the typical size of the region of interest and 
therefore the approximate pathlength of the gamma ray 
through the photon field, riuv is the number density of 
photons near energies corresponding to the peak cross 
section (maximum in a{E, E' ,0)), and a « O.lcrx is this 
approximated cross section. Here ctt is the Thomson 
cross section. 

2.1. Binary companion 

Here we consider the attenuation of gamma rays by 
a high-mass binary companion of the GRB progenitor. 
One class of GRB creation models invokes a compan- 
ion star at a fe w stellar radii to form a rapidly spin- 
ning helium core (Izzard et al.''2004'; 'Podsiadlowski et al.' 
2004: R amircz-Ruiz 2004; Barkov & Komissarov 2010; 
Yoon et al.il20100 . a distance of lO^" '^ to lO"-^ cm for 



a high-mass stars. The effect of radiation from this com- 
panion on gamma-ray flux is limited by the radius R-y 
at which the high energy flux is emitted from the out- 
flow of the progenitor. Most models place this emis- 
sion a t a distance ^ 10^^ cm (Piran' 120041 : IGranot et al.l 
I2008t IKumar fc Barniol Duran .2009i) , typically greater 
than the binary separation. We have then for the opac- 
ity 
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where is the companion luminosity, and is the an- 
gle of photon-photon interactions from Equation ([T]). If 
we consider the opening angle of the jet from relativis- 
tic beaming then the maximum 77 interaction angle 
possible is 

« Rhin/R-i + ^0 Rhin/ Ry + . 

The radius Ry, for models in which the burst is trig- 
gered by the action of the companion, is expected to 
be larger than the progenitor-companion distance i?bin- 
Meanwhile th e high bulk L orentz factor, F ^ 100, in most 
GRB models ()Piranll2004D leads to a tightly beamed radi- 
ation cone. This angular factor therefore greatly reduces 
the likelihood of interaction for GeV-scale GRB emission. 
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Fig. 1. — Minimum gamma-ray production distance from a 
close binary pair at a given separation for which T-y ^ 1. The solid 
and dotted lines arc for gamma-ray energies of 300 and 100 GeV, 
respectively. The lower and upper pairs of curves are for a low- 
metallicity main-sequence companion star of mass 10 Mq and 100 
Mq, see text for details. 

In Figure [1] we show a summary of our results for 
gamma-ray opacity based on binary separation and emis- 
sion radius, assuming that the direction of emission is 
perpendicular to the orbital plane. Lines in this plot 
indicate the minimum gamma-ray production distance 
at which photon-photon opacity is less than unity, for 
a given separation distance between the progenitor and 
its companion. We assume here a low-metallicity main- 
sequence companion at two different stellar masses, 10 
and 1 00 Mf:^, wh ich are modeled as blackbody sources. 
From Pibut et al.l ([1996), we assign these two cases lu- 
minosities of lO'^'^^ and 10^ L0, and effective blackbody 
temperatures of 10^'^ and 10"''^^ K. These results suggest 
that while significant attenuation of gamma-rays can oc- 
cur very close to the progenitor, it is unlikely to be dom- 
inant at the the radius at which gamma-ray are usually 
thought to be produced. 

2.2. Molecular cloud 

For the case of attenuation by the photon field within 
a large molecular cloud, we can look to surveys of these 
complexes within the milky way (MW) and nearby galax- 
ies to estimate the extent and magnitude of a typical 
photon field. Clouds in the MW typically have sizes i?mc 
of lO'^'^ to 1 0^'^ pc, and an H2 gas mass as high as '-^ 10^ 
to 10'^ Mq (jSolomon et al.lll987t iHever et al.il2(30l . and 
other authors confirm similar results for M31 and the 
Large Magellanic Cloud (Rosolowskv 200*3; iFukui et al.l 
I2008f) . With the assumption that 5 to 10% of this gas 
forms stars within a time co mparable to the hfetim e of 
high-mass stars, 20-30 Myrs (|Krumholz et al.l [20061 ). we 
can estimate the highest value likely for our estimate of 
in Equation ^ 



For a Salpeter IMF (|Salpetedll95l 
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we have ~ 0.005 high-mass (> 10 Mq) stars per solar 
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mass of gas converted into stars. If we assume that these 
stars have a typical luminosity of L^, = 10'* Lq which 
is emitted mostly in the UV and is responsible for most 
of this radiation at these wavelengths, then we estimate 
that the most massive molecular clouds could have N^, = 
5 X 10^ massive stars producing N^, = 5 x 10^ Lq in 
optical-UV photons. The energy density in the cloud can 
be approximated as 
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about 6.7 X 10^^ erg cm"'^ under our assumptions. If we 
assume that this energy exists in the form of UV-optical 
(i.e. 3 eV) photons, in a cloud of size R^c = 10 pc, then 
we have an photon number density of 1.4 x 10'^ cm~^. 
Applying Equation ([3]) to find the opacity over a path 
through the cloud, we estimate 

- 3 X 10"^ 

from UV photons in the most massive stellar birth clouds. 
We have not accounted for dust extinction in these sys- 
tems, which can be substantial and would reduce the UV 
flux and t^. 

2.3. Galactic Disk 

While it is generally accepted that long-duration 
GRBs arc associa ted with core-collap se supernovae 
fWooslc v fc Bloomi [20061 iGehrels et all 12009). the re- 
lation between GRB host galaxies and the popu- 
lation of star-forming galaxies appears to be com- 
plex, and recent studies seem to dismiss the idea 
that GR B rates follow star formation i n an unbiased 
manner (iRamirez-Ruiz et all [2002: Cc n fc Fang! 120071 : 
iGuetta fc PiranI 120071 ). Surveys of host galaxies can 
provide a guide for us to choose general parameters in 
this project. A su rvey of GRB host galax ies at a va- 
riety of redshifts bv I Wainwright et al.l ( 20071 ) found that 
hosts can typically be described by an exponential profile 
with a median scale radius of 1.7 kpc; none followed the 
de Vaucouleurs profile associated locally with elliptical 
g alaxies. 

iCastro Ceron et al.l (|2008[ ) analyzed 30 GRB hosts, 
and found typical star-formation rates of 0.01 to 10 
MQ/yr from an assumed unobscured UV, though the in- 
clusion of dust could raise these values considerably. This 
paper also mentions that some galaxies in their sample 
with more poorly constrained SEDs could have star for- 
mation rates as high as a few hun dred MQ/yr. A study 
of 46 hosts bv lSavaglio et al.l (|2009[ ) found star-formation 
rates of 0.01-36 MQ/yr, with a mean of 2.5, after correc- 
tion for dust extinction. Stellar masses in these samples 
ranged from 10^ to 10"'^^ Mq in the former, and 10*'^ to 
1 0^ Mq in t h e latt er. However, 

iPerlev et all (|2009f ) has found that even galaxy hosts 
for obscured "dark" GRBs typically have normal op- 
tical colors overall, suggesting that the dust is un- 
evenly distributed or loca l to the GRB. In their survey, 
ICastro Ceron et al.l (|2008| ) estimat e that ~25% o f GRB 
hosts have dust extinction Ay > 1. iSavaglio et al.l (|2009[ ) 
found an average dust extinction of <Av>~ 0.5 in a sub- 
set of 10 galaxies in their sample; 2 of these galaxies had 
Av > 1. 



Following our method of the last two sections, we 
can perform an order-of-magnitude estimate of for a 
rapidly star-forming galaxy. For a given star-formation 
rate, we can estimate th e UV flux using the approxima- 
tion (jMadau et al.lll998| ) 
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For a star-formation rate of 100 MQ/yr, and a galaxy 
size of 1 kpc, we estimate a total UV output of 2 x 10^^ 
erg/s, and consequently a photon energy density similar 
to that of the molecular cloud in the previous section. 
From Equation ([3]), we estimate a gamma-ray opacity 

- 0.3 

if this energy is emitted as 3 eV photons. 

These finding suggest that the emission from a star- 
forming galaxy, acting over a path length of a kpc or 
more, could be a barrier to GeV gamma-rays from a 
GRB within a galaxy. However, this is only an order- 
of-magnitude estimate of Xy, with a high star-formation 
rate and no account of the absorption and re-emission 
of starlight by dust. Higher star-formation rates than 
we assume here have been seen in GRB-selected sub- 
millimeter galaxies (Micha lowski et al. 2008), and in ex- 
treme cases g alaxies of this type can have rates exceeding 
1000 MQ/yr (jPaddi et al.ll2'005D . albeit with heavy dust 
obscuration. To better understand the circumstances un- 
der which the photon fiux from stars and dust in a galaxy 
might lead to significant attenuation of gamma rays on 
a kiloparsec scale, we have performed a set of numerical 
calculations which are described in the next section. 

3. NUMERICAL CALCULATIONS OF OPACITY FROM 
STELLAR DISKS 

3.1. Model 

In order to accurately estimate the contribution to 
opacity arising from a compact, rapidly star-forming 
galaxy, we have performed a simple simulation of a galac- 
tic disk, and calculated numerically the resulting UV- 
IR photon field and from this the gamma-ray optical 
depth along a given line-of-sight. The stellar spectra 
are calculated usin g the population synthesis code of 
IBruzual fc CharlotJ (120031 ). We assume that the stellar 
population in our galaxy is formed at a constant rate 
over a characteristic timescale Tgf. Stars are then dis- 
tributed in a galactic disk of exponential scale length 
i?e and maximum radius 3i?e, and with thickness 0.3i?e. 
Note that the half-light radius for a exponential profile 
is IMRe. 

Our calculation includes a simple model for dust ab- 
sorption and re-emission . Starl ight is absorbed using the 
model of lCharlot fc Falll ()2000[ ). in which starlight is ab- 
sorbed by 2 components: an interstellar medium (ISM) 
component affecting all stars, and a molecular cloud com- 
ponent affecting only stars younger than 10*" yrs. We 
assume the standard ratio for these opacities, in which 
the ISM is responsible for 30% of the opacity for young 
populations, and the rest arises from absorption within 
molecular clouds. The ISM dust component is assumed 
to have a homogeneous distribution in the same disk as 
the stellar population. All energy absorbed by dust is 
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re-emitted in the IR. We calculate the emiss ion at wave- 
length s above 4 /itm using the templates of iRieke et alj 
(gOOi), which are based on Spitzer observations of local 
galaxies. We consider face-on dust extinction factors of 
0.1, 1.0, and 3.0 in the V-band attenuation ry, broadly 
corresponding to very low, moderate, and large amounts 
of reddening. From this normalization, the extinction 
at other wavelengths is determined by multiplying by a 
power law in wavelength with index -0.7. In galaxies 
with Ty = 3.0, greater than 90% of starlight from stel- 
lar populations younger than 500 Myrs is absorbed by 
dust and re-emitted in the IR, so this corresponds to a 
(U)LIRG-like mode with highly-obscured star- formation. 
All of our galaxy spectra are created using a metallicity 
of [Fe/H] = -1.65. This is somewhat less than the typical 
metallicity of 1 /6 solar that was seen in the sample of 
iSavaglio et al.l (2009); we do not expect this to notably 
impact our final results. 

Figure [5] shows the sample spectra for galaxies of two 
different star- formation timescales (Tgf ) and three values 
for face-on dust extinction (ry = 0.1, 1.0, and 3.0), nor- 
malized to the same mass. For our calculation, we have 
considered timescales over the range 10 Myr < Tgf < 10 
Gyr. In Figure [3l we show a comparison with the SED 
for a starbursting GR B host, GRB000418, as presented 
in iMichalowski et all ((2008) . This source is representa- 
tive of the type of GRB hosts that have the highest in- 
ferred star-formation rate and are therefore expected to 
have the largest opacity to gamma rays. 
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Fig. 2. — The emission spectra for some of the galaxies in our 
model. The solid lines show are for a stellar population with Tgf = 
30 Myr, while the dashed lines are for T^f = 1 Gyr. Line colors 
indicate different levels of dust attenuation, as described in the 
text. Blue, green and red correspond to extinction factors of 0.1, 
1.0, and 3.0 in ry (top to bottom, respectively, at 10^^ Hz). Stellar 
masses have been normalized to 10^^ Mq. 

We do not include a luminosity contribution from AGN 
in our calculation. However, our templates for dust re- 
emission are based on observations of local galaxies, and 
therefore some degree of AGN contamination in the IR 
unavoidably exists. Some degree of evolution in the typ- 
ical IR properties of galaxies likely occurs between local 
redshifts a nd those nea r the p eak of the star-formation 
era (z'^2); IRieke et al.l (|2009( ) estimates that these are 
not larger than a factor of two. In our model we assume 
a scaling law for the galactic disk radius with stellar mass 
of 

Re cx All^^ (8) 
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Fig. 3. — Comparison between one of ou r model spectra and 
that of the fit to GRB 000418 shown in iMTchalowski et al. (200^) 
(black dotted line). Our spectrum (solid red line), is for a total 
dust extinction of ry = 3.0, and has been normalized to match the 
inferred star-formation rate for this host of 288 M0/yr. 

with normalization of 500 pc at 10^ Mq; this slope 
is based upon a rough fit to the half-light radius 
results from a survey of star-forming galaxies in 
iForster Schreiber et al.l ()2009[ ). and is consistent with the 
more compac t galax ies seen in the sample presented in 
iTrujillo et al.l (|2004f ). This particular scaling produces 
galaxies that are more compact than 18 o f 20 GRB hosts 
with rso (w 3i?e) fits in the study of iSvensson et al.l 
((20101) . Since the disk height is also proportional to 
Re in our model, it should be pointed out that for a 
given set of stellar and dust properties the photon flux 
density scales as R~^, and therefore gamma-ray opac- 
ity is proportional to after accounting for the in- 
crease in path length. The UV, optical, and near-IR 
spectral properties of a unit mass of stars are defined in 
our model by the assumed dust attenuation factors and 
star-formation timescale, therefore for gamma-ray ener- 
gies affected solely by these wavelengths ( £ 500 GeV) 
the gamma-ray opacity scales as 

oc Mq/^ (9) 
when these factors are held fixed. 

3.2. Results 

The primary ingredients of our model are the star- 
formation timescale, dust attenuation factor, and galaxy 
size and stellar mass which are conjoined in Equation (jSl), 
as discussed in the previous section. The other impor- 
tant considerations are the location of the GRB within 
the disk and the geometry of the line-of-sight (LOS) path 
from the GRB to the observer. We fix the location of the 
GRB at a radius 0.75i?e (25 % of the simulated radius 
3i?e) and centered within the disk height. We consider 
two possibilities here for the LOS: a "skimming" mode in 
which the path is toward the galactic center at an angle 
of incidence of 10 degrees to the disk, and a geometry in 
which the path is face-on relative to the disk, with angle 
of incidence 90 degrees. The former gives a nearly maxi- 
mal amount of attenuation for a given galaxy, while the 
second gives a closer to average amount of absorption. 

We have calculated opacities for gamma rays from en- 
ergies of 1 GeV to 100 TeV. In Figure [H we show the 
opacity Tj for gamma rays at a particular energy for sev- 
eral sets of galaxy properties. Changing from the skim- 
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ming LOS to one that is perpendicular to the galactic 
disk is found to reduce by a factor of 2 - 3. 
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Fig. 4. — Optical depth for gamma rays created by the photon 
field in and around a galaxy of given stellar properties, for a line-of- 
sight along the galactic disk. Line types and colors are the same as 
in Figure[2] solid and dashed lines are for star-formation timescales 
of 30 Myr and IGyr, and blue, green and red are for dust extinction 
factors of 0.1, 1.0, and 3.0 in tv (top to bottom, respectively, at 
100 GeV). Thin lines for the blue curves show how the attenuation 
changes when the path is instead perpendicular to the galactic 
plane. Results here are again for a galaxy of stellar mass 10^'' Mq. 
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Fig. 5. — The minimum gamma-ray energy for which t-,. > 1, as a 
function of star-formation timescale in the host galaxy. In this plot 
we show results for a low level of dust extinction, Ty = O.l Line 
types indicate various galaxy stellar masses, and we have assumed 
a characteristic radial scale length for the disk from Equation |8] 
Line color indicate results for the different line-of-sight geometries, 
black is for a direction skimming along the disk, and green is for a 
direction perpendicular to the disk, as discussed in the text. The 
thin dotted lines denote star-formation timescales that imply recent 
star-formation rates in excess of 1000 MQ/yr for the given stellar 



The results in Figure [2] for a galaxy of mass 10^° M© 
indicate that attenuation at energies < TeV is signifi- 
cantly less than r-y = 1 even for the optimal "skimming" 
geometry and a fast stellar buildup time of 30 Myrs. As 
mentioned in the previous section, the attenuation is ex- 

2/3 

pected to be approximately proportional to , if we 
assume that galaxy volume scales with mass. The only 
significant attenuation in these cases is due to the popula- 
tion of mid- and far-IR photons created by thermal and 
PAH emission; these being much more numerous than 
the UV, optical, and near-IR photons of direct starlight. 
In Figures [5] and |6] we examine the relationship be- 
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As in Figure [S] but for total dust extinction coefBcient 



tween galaxy mass, star-formation timescale (or star- 
formation rate) and the minimum gamma-ray energy for 
which T.y > 1; these two plots are respectively for cases 
of low and moderate dust extinction (ry — 0.1 and 1.0). 
Results for the largest dust extinction, ry = 3.0, are 
not substantially different from the 1.0 case and are not 
shown. Displaying our results in this way can tell us 
the approximate energy scale at which photon-photon 
opacity becomes an important effect on the high-energy 
spectrum of radiation escaping from the vicinity of the 
galaxy. In general, the opacity only becomes significant 
at multi-TeV energies, where gamma-rays have sufficient 
energy to interact with the re-emitted IR photons. As 
seen in Figure 21 the opacity curve is approximately flat 
between 100 GeV and 1 TeV. Opacity at these energies 
is generally less than 1, and we only find optically-thick 
galaxies in cases of minimal dust extinction, high mass, 
and short star-formation timescale, with star-formation 
rate > 10^ Mo/yr. 

4. DISCUSSION 

Our results suggest that it is unlikely the radiation 
field within a single star-forming region could impact the 
high-energy (GeV and TeV) gamma-ray emission from 
a GRB, although our understanding of these regions is 
limited to nearby galaxies. Significant attenuation due 
to the emission of a massive binary companion is possible 
only if the radius of high-energy photon emission is less 
than predicted by most radiation models. 

Large, rapidly star-forming galaxies can form an 
optically-thick barrier only at multi-TeV energies, where 
attenuation with mid-IR photons can occur, regardless 
of line-of-sight geometry or dust extinction. However, 
photons at these energies would already be strongly at- 
tenuated by extragalactic background light for all but 
the closest GRBs (z ^ 0.05). Creating enough photons 
in the UV, optical and near IR to significantly impact 
transmission of GeV photons is only possible with star- 
formation rates well in excess of 1000 Mq/ji, even with 
a minimal amount of dust attenuation. This requirement 
could be lowered if galaxies are much more compact than 
we have assumed, though the galaxy sizes used here are 
al ready more co n rpact than most of those in the sample 
of iTruiillo et al.l ()2004l ). and these authors mention that 
the densest objects at high redshift tend to be elliptical 
galaxies with much lower photon densities. 

Clumpy or irregular star formation is one possibility 
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that we have not considered in our numerical calcula- 
tions, which assume a smooth exponential disk. The 
total amount of attenuation through the host galaxy is 
proportional to the photon density over the line-of-sight 
path. Therefore the total for the galaxy would only 
be increased significantly from our predictions if a large 
fraction of the total stars in the galaxy were placed in a 
much more compact distribution than we have assumed 
with our exponential profiles. This region would have to 
contain a much higher density of young stars than the 
star- forming regions seen in the MW and LMC, based 
on the order-of-magnitude analysis of §2.21 Merely rear- 
ranging star formation into clumpy regions throughout 
the galaxy would not, in general, affect the total pho- 
ton field or , unless those regions were were positioned 



near the line-of-sight path in a stochastic irregularity. In 
conclusion, both long and short GRB observations can 
be safely used as powerful tools for understanding the 
extragalactic background light . 



We thank J. Granot, X. Prochaska, and J. Primack for 
useful discussions, M. Michalowski and D. Watson for 
providing us with their spectral fits to GRB host galaxies, 
and the anonymous referee for comments which helped 
improve this paper. R.G. acknowledges support from a 
Fermi Guest Investigator Grant during part of this re- 
search. E.R-R acknowledges NSF grant AST-0847563, 
and thanks the Packard Foundation for support. 



REFERENCES 



Abdo, A. A. 2010,'arXiv;1005.0996l 
Abdo, A. A., et al. 2009a, ApJ, 706, L138 
— . 2009b, Science, 323, 1688 

Barkov, M. V., & Komissarov, S. S. 2010, MNRAS, 401, 1644 
Bloom, J. S., Kulkarni, S. R., & Djorgovski, S. G. 2002, AJ, 123, 
1111 

Bruzual, G., & Chariot, S. 2003, MNRAS, 344, 1000 

Castro Ceron, J. M., Michalowski, M. J., Hjorth, J., Malesani, D., 

Gorosabel, J., Watson, D., & Fynbo, J. P. U. 2008, 

arXiv:0803.2235 
Cen, R., & Fang, T. 2007, larXiy:07 10. 43 701 
Chariot, S., & Fall, S. M. 2000, ApJ, 539, 718 
Christensen, L., Hjorth, J., & Gorosabel, J. 2004, A&A, 425, 913 
Daddi, E., et al. 2005, ApJ, 631, L13 
Forster Schreiber, N. M., et al. 2009, ApJ, 706, 1364 
Fruchter, A. S., et al. 2006, Nature, 441, 463 
Fukui, Y., et al. 2008, ApJS, 178, 56 

Gehrels, N., Ramirez-Ruiz, E., & Fox, D. B. 2009, ARA&A, 47, 
567 

Gilmore, R. C, Madau, P., Primack, J. R., Somerville, R. S., & 

Haardt, F. 2009, MNRAS, 399, 1694 
Gilmore, R. C, Prada, F., & Primack, J. 2010, MNRAS, 402, 565 
Gould, R. J., k. Schreder, G. P. 1967, Phys. Rev, 155, 1404 
Granot, J., Cohen- Tanugi, J., & do Couto e Silva, E. 2008, ApJ, 

677, 92 

Guetta, D., & Piran, T. 2007, Journal of Cosmology and 

Astro-Particle Physics, 7, 3 
Heyer, M. H., Carpenter, J. M., k, Snell, R. L. 2001, ApJ, 551, 852 
Iwamoto, K., et al. 1998, Nature, 395, 672 

Izzard, R. G., Ramirez-Ruiz, E., & Tout, C. A. 2004, MNRAS, 
348, 1215 

Krumholz, M. R., Matzner, C. D., & McKee, C. F. 2006, ApJ, 
653, 361 

Kumar, P., & Barniol Duran, R. 2009, MNRAS, 400, L75 
MacFadyen, A. I., Woosley, S. E., & Heger, A. 2001, ApJ, 550, 
410 



Madau, P., & Phinney, E. S. 1996, ApJ, 456, 124 

Madau, P., Pozzetti, L., & Dickinson, M. 1998, ApJ, 498, 106 

Michalowski, M. J., Hjorth, J., Castro Ceron, J. M., & Watson, 

D. 2008, ApJ, 672, 817 
Moskalenko, I. V., Porter, T. A., Strong, A. W. 2006, ApJ, 640, 

L155 

Omodei, N., et al. 2009, GRB Coordinates Network, 9350, 1 
Perley, D. A., et al. 2009, AJ, 138, 1690 
Piran, T. 2004, Reviews of Modern Physics, 76, 1143 
Podsiadlowski, P., Mazzali, P. A., Nomoto, K., Lazzati, D., & 

Cappellaro, E. 2004, ApJ, 607, L17 
Ramirez-Ruiz, E. 2004, MNRAS, 349, L38 

Ramirez-Ruiz, E., Trentham, N., & Blain, A. W. 2002, MNRAS, 
329, 465 

Rieke, G. H., Alonso-Herrero, A., Weiner, B. J., Perez-Gonzalez, 
P. G., Blaylock, M., Donley, J. L., & Marcillac, D. 2009, ApJ, 
692, 556 

Rosolowsky, E. 2007, ApJ, 654, 240 
Salpeter, E. E. 1955, ApJ, 121, 161 

Savaglio, S., Glazebrook, K., & Le Borgne, D. 2009, ApJ, 691, 182 
Solomon, P. M., Rivolo, A. R., Barrett, J., & Yahil, A. 1987, ApJ, 
319, 730 

Stecker, F. W., de Jager, O. C, Salamon, M. H. 1992, ApJ, 390, 
L49 

Svensson, K. M., Levan, A. J., Tanvir, N. R., Fruchter, A. S., 

Strolger, L. G. 2010, MNRAS, 405, 57 
Tout, C. A., Pols, O. R., Eggleton, P. P., & Han, Z. 1996, 

MNRAS, 281, 257 
Trujillo, I., et al. 2004, ApJ, 604, 521 

Wainwright, C, Berger, E., & Penprase, B. E. 2007, ApJ, 657, 367 
Woosley, S. E. 1993, ApJ, 405, 273 

Woosley, S. E., & Bloom, J. S. 2006, ARAfc A, 44, 507 

Yoon, S., Woosley, S. E., & Langer, N. 2010. larXiv:1004.0843l 



